ADA  0  8  66  7  1 


A  REVIEW  OF  NONLINEAR  FRACTURE 
MECHANICS  RELATIVE  TO  FATIGUE 


T.  Weerasooriya 
J.  P.  Gallagher 
H. .  C .  Rhee 


University  of  Dayton  Research  Institute 
300  College  Park  Drive 
Dayton ,  Ohio  45469 


S£.pT£/*B£X  3°,  /97? 


TECHNICAL  REPORT  AFML-TR-79-4196 

Final  Report  for  Period  September  1978  -  September  1979 


Approved  for  public  release;  distribution  unlimited 


AIR  FORCE  MATERIALS  LABORATORY 
AIR  FORCE  WRIGHT  AERONAUTICAL  LABORATORIES 
AIR  FORCE  SYSTEMS  COMMAND 
WRIGHT-PATTERSON  AFB,  OHIO  45433 


DTIC 


SELECTE 
JUL  1  1  1980 


B 


NOTICE 


When  Government  drawings,  specifications,  or  other  daia  are  used  for  any 
purpose  other  than  in  connection  with  a  definitely  related  Government 
procurement  operation  the  United  States  Government  thereby  incurs  no 
responsibility  nor  any  obligation  whatsoever;  and  the  fact  that  the  government 
may  have  formulated,  furnished,  or  in  any  way  supplied  the  said  drawings, 
specifications,  or  other  data,  is  not  to  be  regarded  by  Implication  or 
otherwise  as  in  any  manner  licensing  the  holder  or  any  other  person  or 
corporation,  or  conveying  any  rights  or  permission  to  manufacture,  use,  or 
sell  any  patented  invention  that  may  in  any  way  be  related  thereto. 

This  report  has  been  reviewed  by  the  Information  Office  (01)  and  is 
releasable  to  the  National  Technical  Information  Service  (NTIS).  At  NTIS,  it 
will  be  available  to  the  general  public,  including  foreign  nations. 

This  technical  report  has  been  reviewed  and  is  approved  for  publication. 


THEODORE  NICHOLAS 
Project  Engineer 


If}  _ 

NATHAN  G.  H’I/peR,  Chief 
Metals  Behavior  Branch 
Metals  and  Ceramics  Division 


If  your  address  has  changed,  if  you  wish  to  be  removed  from  our  mailing  list, 
or  if  the  addressee  is  no  longer  employed  by  your  organization  please  notify 

_ ,  W-PAFB,  OH  45433  to  help  us  maintain  a  current  mailing 

list. 


Copies  of  this  report  should  not  be  returned  unless  return  is  required  by  se¬ 
curity  considerations,  contractual  obligations,  or  notice  of  a  specific  document 

AIR  FORCE/88780/17  Juna  X980  -  800 


UMf:T,ASfiTF'i-E:n 


SECURITY  CLASS'FIC.ATION  OF  THIS  “tss  I'tFh-fi  D*t»  Enlrrrd) 


■/ 


M 


Q  ■? JREPORT  DOCUMENTATION  PAGE 

.  J*1  macn - - - - r rrr-r 


I  REPORTtKlMQEn 


'Qj 


H 


l  AFMIMTR-79-4196 ; 

<"  T I  T_L 

A  Review  of  Nonlinear  Fracture 
Mechanics  Relative  to  Fatigue, 


!■>  GOVT  ACCESSION  NO. 


RFFP  INSTRUCTIONS 
Otf  ORT  COMPLETING  FORM 


»•  ^E&cj?:x±sit>XAf  nr,  number 

■  •  Pi 


7 


«Keari>iWW«WMCovER 

a«w«rL  RepSSt  ,  T - 

Sep#  78— ."Sep#  79„ 


7.  .avthor£»,>_ 


/'  *  j  ,  j  y-  !.f 

— - T  Ut  Jr  udr-tr-79-9/1/ 

-  ■_■  t — — number/s; 


T.  I Weerasooriya ,  if.P./  Gallagher/'  H.C.j Rhee^  '  F 33615-7 8-C-5184 

f..  .  ■/ — “K— — — ,  \j£T’ "  "  -—v-—.—.  •— /- 

t^5pr\* 


9.  PERFORMING  0R(  ANlZATlON  NAME  AND  ADDRESS 

university  of  Dayton  V  ;v 
Research  lnsta*»twi  ^ 

Dayton ,  Ohio  45469  _ 


'  J 


M  '  ..’..MEMT.  PROJECT,  TASK 
*  R-'  INIT  NUMBERS 


II.  CONTROLLING  OFFICE  NAME  AMO  ADDRESS 

Air  Force  Materials  Laboratory 
Air  Force  Systems  Command 
Wright-Patterson  AFB ,  Ohio  45433 


am 


U.  MONITORING  AGENCY  NAME  6  ADDRESS^/, dl/(e,onf  from  Controlling  Ollic*T 


Sepl 

r.frr^'gfrrgyinrr 


lzE7 


76 


15.  SECURITY  CLASS,  (of  thi,  wpo.v) 

Unclassified 


15a.  DECLASSIFICATION  DOWNGRADING 
SCHEDULE 


16.  ‘distribution  statement  fof  thU  R»pt.rt) 

Approved  for  public  release;  distribution  unlimited, 


17.  DISTRIBUTION  STATEMENT  (of  the  «b»fra-f  en farad  In  Block  JO,  it  different  from  Report) 


IS.  SUPPLEMENTARY  NOTES 


19*  KEY  WOROS  fContfmi*  on  tevetee  eide  if  neceeeery  end  idbntlfy  by  block  number) 


20.  ABSTRACT  (Continue  on  revree  elde  If  neceeemry  end  identify  by  block  number) 


t  I 


The' report  reviews  currently  available  parameters,  methodology. 
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fracture  mechanics  (EPFM)  technology  for  modeling  fatigue  crack 
growth  rate  (FCGR)  behavior.  In  general,  it  can  be  stated  that 
linearly  elastic  fracture  mechanics  (LEFM)  parameters  can  be 
used  to  model  FCGR  behavior  except  when  the  zone  of  cyclic 
plasticity  accompanying  the  crack  process  is  large  with  respect 
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SECTION  I 


INTRODUCTION 

1.1  BACKGROUND  AND  SCOPE 

In  the  mid  1950's,  a  new  concept  for  characterizing  the 
fracture  behavior  of  materials  evolved  from  studies  of  the  stress 
field  at  the  tip  of  a  crack  in  an  elastic  body.  It  was  quickly 
noted  that  a  single  parameter,  the  stress  intensity  factor  (K) , 
incorporated  the  influence  of  remote  stress,  crack  length,  and 
structural  geometry.  This  parameter  was  successfully  used  to 
correlate  the  conditions  associated  with  the  initiation  of  fracture 
between  service  structures  and  laboratory  test  specimens  manufactured 
from  high  strength  metals.  By  the  early  1960's,  this  same  concept 
was  found  to  provide  an  acceptable  procedure  for  correlating 
fatigue  crack  growth  rate  behavior  between  cracks  of  different 
configurations  subjected  to  various  loading  conditions.  The 
approach  to  describing  fracture  and  fatigue  crack  growth  rate 
behavior  based  on  the  stress  intensity  factor  has  sometimes  been 
referred  to  as  the  fracture  mechanics  approach,  but  more  specifi¬ 
cally  and  recently  as  the  linear  elastic  fracture  mechanics  (LEFM) 
approach. 

The  LEFM  approach  basically  follows  a  similitude  approach 
where  the  identical  local  stress-strain  behavior  and  resulting 
crack  growth  behavior  are  assumed  to  occur  for  all  cracks  having 
the  same  stress  intensity  factor.  It  is  realized  that  ordinary 
structural  materials  do  not  exhibit  complete ' elastic  behavior, 
especially  in  the  crack  tip  region.  However,  if  the  zone  of 
nonlinear  deformation  is  small  and  contained  in  the  elastic  stress 
field,  the  elastic  stress  field  will  set  the  boundary  conditions 
for  the  deformation  of  the  material  in  the  nonlinear  zone.  Hence, 
the  LEFM  approach  will  describe  fracture  and  fatigue  crack  growth 
rate  behavior  as  long  as  the  nonlinear  deformation  zone  at  the 
crack  tip  is  small  enough  so  that  the  surrounding  elastic  stress 
field  can  be  described  by  the  elastic  stress  intensity  factor. 


Since  1970,  an  ever  increasing  amount  of  attention  has  been 
given  to  defining  the  limitations  of  the  widely  used  LEFM  approach. 
Most  of  this  attention  is  devoted  to  improving  the  characterization 
of  the  onset  of  monotonic  cracking  behavior  and  the  subsequent 
slow  stable  tearing  that  precedes  the  fracture  of  ductile 
materials  (of  the  type  that  might  be  found  in  nuclear  power  plant 
systems) .  More  recently,  loading  and/or  material  conditions  have 
been  noted  to  limit  the  LEFM  parameter  K's  ability  to  correlate  the 
fatigue  crack  growth  rates  generated  using  (a)  different  types  of 
laboratory  specimen  geometries  and  (b)  different  stress-crack 
length  combinations  within  the  same  specimen  geometry.  Unfortunately, 
the  conditions  that  cause  the  normally  good  fatigue  crack  growth 
rate  (FCGR)  correlations  (based  on  the  LEFM  approach)  to  breakdown 
are  not  particularly  well  documented.  Typically,  when  the  FCGR- 
LEFM  correlations  are  noted  to  breakdown,  the  resulting  behavior 
is  referred  to  as  nonlinear  fracture  mechanics  behavior.  This 
terminology,  however,  has  led  to  even  further  communication 
problems . 

This  report  has  been  written  in  an  attempt  to  clarify  the 
limitations  of  the  FCGR-LEFM  methodology  and  to  review  alternate 
approaches  which  have  the  potential  for  improving  the  characteriza¬ 
tion  of  FCGR  behavior.  Additionally,  the  report  provides:  (1)  a 
review  of  currently  available  elastic-plastic  fracture  mechanics 
(EPFM)  parameters  that  might  improve  FCGR  correlations,  (2)  a 
summary  of  available  data  that  exhibit  a  lack  of  correlation  based 
on  the  LEFM  approach,  and  (3)  a  recommended  course  of  action  for 
investigating  parameters  that  have  the  potential  to  provide  for 
improved  FCGR  correlations. 

1.2  THE  BREAKDOWN  OF  LEFM 

1.2.1  The  LEFM  Approach  to  Fatigue 

To  clarify  the  limitations  in  the  LEFM  approach,  we 
will  initially  provide  an  overview-type  review  of  the  FCGR-LEFM 
methodology.  Shown  in  Figure  1  is  a  schematic  which  describes 
the  procedure  by  which  fatigue  crack  growth  (FCG)  data  are 
collected  in  the  laboratory,  reduced  to  FCGR  data,  and  subsequently 
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Figure  1.  Schematic  Illustrating  FCGR  Data  Reduction 
Approach  Based  on  the  LEFM  Parameter,  K 
The  Maximum  Stress  Intensity  in  a  Stressax 
Cycle  (Reference  1) . 
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portrayed  as  a  function  of  the  LEFM  parameter  K.  To  summarize  the 
approach,  FCG  data  are  collected  under  constant  amplitude  stress 
cycling  histories  at  various  stress  ratios  using  a  well  characterized 
structural  geometry,  i. .e.  characterized  with  respect  to  the  crack 
driving  parameter.  The  FCG  data  are  reduced  using  standardized 
numerical  techniques  and  presented  as  a  function  of  the  LEFM  crack 
driving  parameter . 

Various  authors  choose  to  use  K  ,  the  maximum  value 

max 

of  the  stress  intensity  factor  (SIF)  in  a  cycle  of  loading,  or  AK, 
the  SIF  range  (maximum  minus  minimum  SIF)  as  the  LEFM  parameter. 

ASTM  Subcommittee  E24,04  on  "Subcritical  Crack  Growth"  has  adapted 
AK  as  the  parameter  that  will  be  used  in  portraying  FCGR  data 
generated  according  to  the  new  ASTM  Standard  E647-78T. 

The  usefulness  of  the  FCGR-LEFM  data  comes  from  the 
observation  that  these  data  are  independent  of  both  structural 
geometry  (including  crack  size)  and  stress  level..  In  fact,  this 
independence  is  the  basic  hypothesis  underlying  what  is  referred  to 
as  the  fracture  mechanics  approach.  Numerous  instances,  where  the 
FCGR-LEFM  modeling  approach  has  been  suggested  to  fail ,  have 
resulted  from  less  than  accurate  characterizations  of  the  SIF; 
this  is  particularly  the  case  when  one  attempts  to  model  3-D  flaws 
in  stress  concentration  regions  with  gross  estimates  of  the  SIF. 

As  suggested  by  Figure  2,  the  process  of  numerical 
differentiation  of  FCG  data  is  reversed  when  one  desires  to  obtain 
a  life  prediction  for  a  structural  component  which  contains  a  crack. 
The  final  verification  for  any  life  prediction  scheme  comes  when 
life  predictions  are  made  for  conditions  (stress  level,  stress  ratio, 
structural  geometry,  etc.)  for  which  the  test  life  is  not  known 
a  priori  to  the  analyst;  this  is  normally  referred  to  as  making  a 
"blind"  life  prediction.  In  general,  FCG  blind  life  predictions 
within  plus  or  minus  twenty  (20)  percent  of  the  observed  test  life 
are  considered  good  while  those  that  differ  by  more  than  a  factor 
of  two  (2)  are  considered  poor. 
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1.2.2  High-FCGR  Nonlinear  Behavior 


Assuming  that  the  LF.FM  approach  applies,  material 
engineers  can  expect  the  FCGR  data  generated  using  any  well 
characterized  (cracked)  structural  geometry  to  always  scatter 
closely  about  the  mean  FCGR  behavior  curve  established  from  data 
collected  using  any  other  well  characterized  structural  geometry. 

One  of  the  initial  definitions  of  what  might  be  called  nonlinear 
fracture  mechanics  behavior  arose  as  a  result  of  portraying  FCGR 
data  scattered  around  a  power  law  equation,  attributed  to  Paris2, 
and  having  the  form: 

gf=C(AK)n  (1) 

which  describes  a  straight  line  in  log-log  coordinates.  A  set 
of  power  law  equations,  in  terms  of  Kmax  can  be  seen  expressed  in 
Figures  1  and  2.  Figure  3  provides  a  complete  representation  of 
FCGR  data  for  a  single  stress  ratio  (R)  that  fully  illustrates 
the  shape  that  a  FCGR-LEFM  curve  assumes.  As  the  reader  will  note, 
a  power  law  equation  will  have  its  greatest  utility  in  the  central 

3 

region  of  the  FCGR-LEFM  curve.  Barsom  ,  in  considering  the  strain 
energy  release  rate  (AG  =  AK2/E)  approach  to  fatigue,  observed 
that  deviations  from  the  power  lav/  relation,  which  he  termed 
"nonlinear  behavior",  could  be  anticipated  on  the  basis  of  a  crack 
opening  displacement  (COD)  parameter.  Barsom  used  LEFM  methodology 
to  estimate  the  controlling  COD  parameter  from  the  S1F  and  was  not 
strictly  interested  in  the  limitations  of  LEFM. 

Additional  investigators  who  have  also  studied  FCGR 
behavior  above  the  "power  law"  region  of  the  FCGR-LEFM  curve  have 
noted  that  FCGR  behavior  is  not  always  uniquel}  related  to  the  LEFM 

4 

driving  parameter.  Frost  et  al.  have  presented  previously  published 
FCGR  data  as  a  function  of  AK  which  demonstrate  that  the  maximum 
applied  stress  (or  net  section  stress)  will  influence  the  FCGR 
behavior  in  the  High-FCGR  region  for  some  but  not  all  metals 
studied.  Dubensky^  ,  in  his  analysis  of  7075  and  2024  aluminum 
alloys  in  the  High-FCGR  region, has  noted  a  controlling  effect  of 
the  maximum  applied  stress  on  the  2024  behavior  (See  Figure  4) 
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Figure  3.  FCGR-LEFM  Curve  for  Complete  Range  of 
Observed  Behavior. 
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Figure  4 


LEFM-AK  Parameter  Correlation  with  FCGR  Data 
for  2024-T3  Aluminum  Alloy  (Reference  5) . 


but  not  on  the  7075.  Dowling  and  Walker**  have  presented  data  for 
two  relatively  low  strength  steels,  Man-Ten  and  RQC-100.  They 
have  noted  evidence  for  stress  or  crack  length  effects  only  in  the 
High-FCGR  region  in  the  lower  strength  Man-Ten  steel  as  shown  in 
Figure  5.  While  at  this  time  it  can  not  be  stated  absolutely, 
it  does  appear  that  the  FCGR  behavior  of  the  more  ductile  alloys 
exhibit  an  effect  of  maximum  stress  in  the  High-FCGR  region,  thus 
violating  a  basic  premise  of  the  LEFM  approach  to  fatigue.  For 
the  brittle  alloys,  this  effect  may  not  be  observed  because  fracture 
takes  place  before  the  plastic-zone  can  grow  to  a  size  that 
invalidates  the  LEFM  assumptions  of  contained  plasticity.  Frost 
et  al.^  attempt  to  qualify  the  usefulness  of  LEFM  by  setting  a 
limit  on  the  plastic  zone  size  (r^)  in  a  center-cracked  (CC)  or 
edge-cracked  (EC)  test  specimen: 


r 


y 


V 

max  ^  a 

a  '  -  7 

o 


(2) 


where  oq  is  the  yield  strength  and  a  is  the  half-crack  length  in 
the  CC  specimen  and  the  crack  length  in  the  EC  specimen.  While 


Frost  etal.  support  their  hypothesis  with  High-FCGR  observations 
for  a  number  of  materials,  their  testing  programs  were  not  designed 
to  validate  the  use  of  Equation  2  as  an  absolute  condition  for 


applying  the  LEFM  approach. 


1.2.3  Nonlinear  Behavior  at  Notches 


The  material  at  the  root  of  a  notch  in  a  highly 
loaded  structure  does  not  need  to  contain  a  pre-existing  defect 
to  experience  fatigue  cracking  problems.  When  the  level  of  cyclic 
plastic  deformation  is  consistently  high,  a  crack  will  initiate 
and  subsequently  grow  under  the  influence  of  this  deformation. 
Hammouda  and  Miller7  have  recently  provided  a  clear  description 
of  the  expected  FCGR  behavior  in  notch  regions  when  the  level  of 
cyclic  loading  induces  relatively  large  local  plastic  deformation. 
They  accurately  point  out  that  LEFM  technology  is  invalid  here 
because  the  crack  is  imbedded  within  a  plastic  stress-strain  field 
and  not  an  elastic  stress  field  generated  by  the  moving  crack. 
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Figure  5.  LEFM-AK  Parameter  Correlations  with  FCGR  Data 
for  Man-Ten  Steel  (Reference  6) ;  The  Parameter 
AK  is  Defined  as  Kmax  -  Kminr  where  Kmin  =  amin-  (K./o) 
The  Data  Were  Generated  Using  WOL  Specimens. 


Figure  6  portrays  the  FCGR  pattern  observed  by  the  small  cracks 

initiated  by  and  embedded  within  the  cyclic  plastic  stress-strain 
/ 

field  at  the  root  of  a  notch.  As  the  figure  shows,  a  crack  first 
grows  at  rapid  rates  and  as  its  length  increases , the  rate  of  growth 
decreases  until  the  behavior  becomes  controlled  by  the  assumptions 
of  LEFM.  Based  on  a  modification  suggested  by  El  Haddad  et  al.  , 
it  may  be  possible  to  distinguish  between  cracks  that  can  be 
modeled  using  a  stress-strain  field  parameter  such  as  the  J-Integral 
and  those  that  might  require  parameters  that  are  beyond  those 
currently  available. 

1.3  SUMMARY  OF  NONLINEAR  BEHAVIOR 

The  correlation  of  FCGR  with  the  single  LEFM  parameter 
appears  to  fail  for  those  conditions  in  which  the  crack  is  either 
propagating  into  or  through  a  zone  of  gross  plastic  yielding 
(Figures  4-6) .  When  a  crack  is  propagating  through  a  material 
experiencing  net  section  yielding,  material  ratcheting  will  further 
increase  the  levels  of  crack  tip  strain  fields  above  that  pre¬ 
dicted  by  LEFM.  In  these  cases,  the  crack  no  longer  generates 
an  elastic  stress  field  which  contains  and  controls  the  cyclic 
behavior  of  the  nonlinear  material  at  the  tip  of  the  crack. 

To  better  describe  FCGR  response  under  loading  conditions 
which  induce  nonlinear  material  behavior,  a  number  of  investigators 
have  suggested  the  use  of  Elastic-Plastic  Fracture  Mechanics 
(EPFM)  parameters.  These  parameters  are  based  on  either  the 
field  approach  or  the  crack  tip  approach.  The  field  approach 
follows  that  of  the  LEFM  approach  in  that  a  parameter  is  used  that 
can  be  related  to  the  intensity  or  magnitude  of  the  response-control¬ 
ling  elastic-plastic  stress-strain  field  at  the  tip  of  a  crack.  The 
crack  tip  approach  is  based  on  a  philosophy,  whereby  one  determines 
the  value  of  the  parameter  at  the  crack  tip  which  is  suspected 
of  controlling  the  crack/growth  behavior.  When  EPFM  parameters 
are  used  to  describe  the  fatigue  and  fatigue  crack  growth  behavior, 
the  approach  is  normally  referred  to  as  the  nonlinear  fracture 
mechanics  (NLFM)  or  elastic-plastic  fracture  mechanics  (EPFM) 
approach;  its  characterizing  parameter  is  referred  to  as  the  NLFM 
or  EPFM  parameter. 
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TOTAL 


Figure  6 


Plastic  and  Elastic  Fracture  Mechanics 
Characterization  of  Fatigue  Crack  Growth 
(Reference  7) . 


In  the  sections  that  follow,  we  will  review  currently 

available  EPFM  parameters  and  FCGR  data  which  support  the 

development  of  a  FCGR-EPFM  metnodology  along  the  lines  suggested 

by  Figures  1  and  2  where  K  „  is  replaced  with  an  EPFM  parameter  "P". 

max 

The  basic  restriction  of  such  a  development  is  that  the  FCGR- 
EPFM  curve,  once  established,  will  be  independent  of  both  the  level 
of  applied  cyclic  stress  and  structural  geometry  (including  crack 
length) . 
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SECTION  II 

NONLINEAR  FRACTURE  MECHANICS  PARAMETERS 

A  survey  of  potentially  useful  parameters  applicable  for 
FCGR  correlations  has  shown  that  the  available  parameters  can  be 
divided  into  two  categories: 

(1)  those  that  characterize  the  intensity  of  an  elastic- 
plastic  stress-strain  field,  and 

(2)  those  that  characterize  the  magnitude  of  a  suspected 
controlling  parameter  at  a  point  located  at  or  near  the 
crack  tip. 

This  section  identifies  the  parameters  by  category  and  provides 
the  background  associated  with  the  development  of  these  parameters. 
Most  of  the  parameters  identified  in  this  section  have  only  been 
used  so  far  in  the  monotonic  fracture  loading  situations.  Also, 
currently  available  numerical  techniques  for  evaluating  some  of 
these  parameters  are  reviewed. 

2.1  THE  FIELD  APPROACH 

There  are  a  broad  class  of  parameters  which  have  evolved 
from  analytical  studies  of  the  crack  tip  stress  and  strain  fields. 
The  analytical  continuum  mechanics  studies  have  shown  that  single 
parameter  characterizations  of  the  stress-strain  field  at  a  crack 
tip  are  possible  as  long  as  these  fields  are  not  distorted  by 
gross  plastic  deformation.  It  is  recognized  that  these  stress- 
strain  fields  do  not  extend  directly  into  the  intensely  distorted 
localized  region  at  the  crack  tip,  where  the  deformation  processes 
are  controlling  the  observed  crack  growth  behavior.  However, 
since  the  intensely  distorted  localized  region  is  imbedded  within 
a  less  intense  but  definable  stress-strain  field,  it  is  argued 
that  the  less  intense  stress-strain  field  provides  the  boundary 
conditions  that  control  the  deformation  processes  within  the 
intensely  distorted  zone.  Thus,  if  the  parameter  which  establishes 
the  magnitude  of  the  less  intense  stress-strain  field  can  be 
determined,  it  can  be  utilized  via  a  similitude  hypothesis  to 
characterize  the  behavior  of  cracks. 
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This  subsection  will  identify  the  currently  available  field 
parameters.  Table  1  provides  a  listing  of  those  parameters  which 
have  been  used  to  characterize  the  intensity  of  the  crack  tip 
stress-strain  field.  Several  LEFM  derivatives  are  presented 
in  the  table  for  completeness. 


TABLE  1 

NONLINEAR  FRACTURE  MECHANICS  FIELD  PARAMETERS 


Elastic  Field  Derivatives 

Stress  Intensity  Factor  (SIF) 
Plasticity  Corrected  SIF  (a*=a+ry) 
Strain  Intensity  Factor 
Neuber  Intensity  Factor 

Elastic-Plastic  Fields 

J-Integral 

Crack  Tip  Opening  Displacement  (CTOD) 


2.1.1  Derivatives  of  Elastic  Fields 

With  the  success  of  LEFM  in  predicting  extension 
of  cracks  in  elastic  materials,  researchers  have  been  considering 
the  extension  of  LEFM  to  the  elastic-plastic  regime. 


Figure  7  shows  the  elastic-stress-field  distribution 
ahead  of  the  crack  tip  when  there  is  no  plasticity  (a  case 

which  does  not  exist  for  real  materials)  the  Mode  I  crack  opening 
stress  oy  is  proportional  to  the  stress  intensity  factor  (SIF),  Kj, 
in  the  vicinity  of  the  crack  tip.  The  stress-field  equations 
show  that  the  distribution  of  the  elastic-stress  field  in  the 
vicinity  of  the  crack  tip  is  invariant  in  all  structural  components 
subjected  to  this  type  of  deformation,  and  that  the  single 
parameter  Kj  describes  the  intensity  of  crack  tip  stress  and  strain 
field.  The  rationale  in  using  the  stress  intensity  factor  to 
describe  FCGR  behavior  follows  from  the  ability  of  the  elastic 
stress-strain  field  to  control  the  deformation  behavior  in  the 
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Figure  7 


Elastic-Stress-Field  Distribution  Ahead  of 
Crack  (Reference  9) . 


distorted  crack  tip  zone  under  cyclic  loading.  This  same 
rationale  holds  for  the  strain  intensity  factor  (K£). 

In  metals,  stresses  at  the  crack  tip  are  finite  and 
plastic  yielding  occurs  at  the  crack  tip.  A  schematic  illustrating 
the  distribution  of  the  stress  normal  to  the  path  of  the  crack  is 
shown  in  Figure  8.  The  extent  of  the  yielding  ahead  of  the  crack 
tip  can  be  estimated  as  twice  the  Irwin  plastic  zone  size  (r^; 
which  is  given  by: 


for  plane-stress  conditions  and 


1  KI  2 

V  (3b» 

for  plane-strain  conditions.  Irwin^-2  suggested  that  an  increase 
in  the  real  crack  length  ry  (See  Figure  9) 13  will  result  in  a 
modified  crack  which  has  an  elastic  stress  distribution  similar 
to  the  elastic- field-stress  distribution  in  the  vicinity  of  the 
real  crack's  plastic  zone  boundary.  Therefore,  when  the  amount 
of  yielding  at  the  crack  tip  is  limited ,  attempts  have  been  made 
to  use  LEFM  parameters,  corrected  to  account  for  the  plasticity, 
to  modeling  FCGR  behavior.  In  addition  to  the  plastic  zone 
correction  suggested  by  Irwin,  there  exist  other  corrections  by 
Dugdale  H  and  Newman  15. 

Attempts  to  model  FCGR  behavior  using  plasticity 
modified  SIF  for  elastic-plastic  or  fully-plastic  material  behavior 
have  had  only  modest  success.  The  reason  for  the  limited  success 
results  from  the  fact  that  when  the  plastic  zone  size  is 
appreciably  large,  (larger  than  the  K-field;  zone  1  in  Figure  9) 
the  stress-field  in  the  vicinity  of  the  crack  tip  is  not  predicted 
accurately  by  the  elastic  solution.  Thus,  for  more  generality 
and  accuracy  one  must  proceed  to  an  elastic-plastic  analysis; 
that  is, to  consider  the  vicinity  of  the  crack  tip  as  an  elastic- 
plastic  or  fully  plastic  field. 


Figure  9.  Linear-Elastic  Crack-Tip  Stress  and  Strain 
Fields  (Reference  13) . 

2.1.2  Elastic-Plastic  Field  Parameters 

The  solutions  for  the  in-plane  tensile  opening  mode 
(Mode  I  type)  of  deformation  problems  are  the  primary  interest 
in  the  present  work  as  well  as  fracture  mechanics  in  general. 
However,  the  mathematical  difficulties  have  prevented  detailed 
treatment  of  elastic-plastic  problems.  Except  in  some  out-of- 
plane,  tearing  mode  (Mode  III  type)  of  problems1^ ' ^ ,  rigorous 
mathematical  solutions  of  elastic-plastic  problems  are  not 
available  in  general;  the  available  limited  cases  of  elastic- 
plastic  crack  tip  stress  analyses  will  be  reviewed  in  this 
paragraph. 

Figure  10  shows  the  crack  tip  and  area  ahead  of 

]  O  , 

the  crack  tip  .  The  region  ahead  of  the  crack-tip  is  divided 
into  three  distinct  zones:  (1)  elastic,  (2)  elastic-plastic, 
and  (3)  intensely  nonlinear  (large  strains  and  rotations, 
and  ductile  cavities)  zone.  The  elastic  zone  (1)  controls  the 
behavior  when  the  plastic  zone  size  is  small  compared  to  the 
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Zone  1  =  An  Elastic  Field  Surrounding  the  Crack  Tip 

Zone  2  =  An  Elastic-Plastic  Field  Surrounding  the  Crack  Tip 

Zone  3  =  An  Intense  Zone  of  Deformation 

Figure  10.  Crack-Tip  Stress  and  Strain  Fields 

Surrounding  the  Crack  (Reference  13) . 

elastic  boundaries.  In  this  case,  referred  to  as  small  scale 
yielding,  LEFM  is  applicable.  If  the  plastic  zone  size  is  large, 
compared  to  the  case  of  small  scale  yielding,  LEFM  is  not  applicable. 

The  intense  elastic  plastic  stress-strain  field 
contained  within  zone  2  of  Figure  10  is  further  expanded  in 
Figure  11. 


Figure  11.  Expanded  View  of  the  Elastic-Plastic 
Stress-Strain  Field  (Reference  13) . 
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When  the  intensely  deformed  process  zone  is  small  compared  to 

the  size  of  the  elastic-plastic  zone  under  consideration,  the 

deformation  theory  of  plasticity  for  a  power  hardening  material 

can  be  used  to  obtain  stress-strain  solutions  ahead  of  the  crack 

tip  outside  the  intensely  deformed  process  zone  as  suggested  by 
18  19 

Hutchinson  and  Rice  and  Rosengren  .  These  authors  expressed 
power  hardening  using  a  stress  (a) -strain  (e)  relationship  that 
is  given  by: 

'  a  -  a  (-§-  )N  (4) 

o 


where  oQ  and  sq  are  the  flow  stress  (yield  strength)  and  strain, 
respectively,  and  N  is  the  strain  hardening  exponent. 

When  Equation  4  is  used  to  model  the  behavior  of 
the  material  in  the  plastic  range,  the  stress,  strain,  and 
displacement  functions  for  the  crack  tip  region  are  given  by 
Equation  5,  regardless  of  the  amount  of  plastic  deformation: 
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(5a) 

(5b) 


1 

=  Kr  •  u^QfN)*  r  N+1  (5c) 

where  and  K£  are  the  magnitudes  of  singularities  of  appropriate 

quantities  with  K£  being  a  function  of  .  The  functions  5^, 

e. .  and  CL  depend  on  angle  and  exponent  N  in  Equation  4.  Equations 

1  j  1 

5a, 5b  and  5c  has  been  referred  to  as  the  set  of  "HRR"  field  equations 
after  the  initial  investigators  18,19. 

In  the  derivation  of  the  HRR  field  equations,  the 
deformation  theory  of  plasticity  was  used  to  describe  the  material 
behavior.  Thus,  when  unloading  occurs  during  deformation,  the 
loading  path  is  different  from  the  proportional  loading  assumed 
and  the  validity  of  the  above  solutions  is  not  guaranteed. 
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The  magnitude  of  the  singularity  in  Equation  5 
in  terms  of  J  can  be  written  as^ 


(6a) 

oij  (0,N) 

(0,N)  (bb) 

N 

rN+1  ui(0,N)  (6c) 


when  the  process  zone  is  small.  The  J  parameter  is  given  by  the 
contour  integral20 

J  =  wdy  -  *•  ds  (7) 

Jr 

which  has  been  determined  to  be  independent  of  the  location  of 
the  contour  curve  r . 

Hence,  for  a  given  material,  with  the  assumption  of 
power  law  hardening  behavior,  deformation  theory  of  plasticity, 
and  proportional  loading,  there  exists  a  unique  elastic-plastic 
stress  and  strain  field  which  is  characterized  by  its  intensity  J. 

In  a  manner  similar  to  the  LEFM  approach  where  K  measures  the  intensity 
of  stress  and  strain  within  the  elastic  crack  tip  field,  the  para¬ 
meter  J  defines  the  intensity  of  the  elastic-plastic  stress  and 
strain  at  the  crack  tip  field  and  thus  provides  a  basis  for  a 
nonlinear  fracture  mechanics  approach.  The  use  of  J  to  define  the 
level  of  elastic-plastic  stresses  and  strains  around  the  crack  tip 
requires  that  the  intensely  deformed  process  zone  is  small. 
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For  ideally  plastic  materials,  the  relationship 

between  the  J-Integral  and  the  crack  tip  opening  displacement 

(defined  by  the  opening  distance  between  the  intercepts  of  two  45° 

lines,  drawn  back  from  the  crack  tip,  with  the  deformed  profile  of 

20  21 

the  stationary  crack)  is  given  by  ' 


J  =  a  aQ  6fc  (8) 
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where  a  is  a  constant.  More  recently  Shih  and  co-workers  ' 
have  shown  that  Equation  8  can  be  used  to  relate  J  and  for 
strain  hardening  materials  where  the  constant  a  is  replaced  with 
a  function  which  is  strongly  dependent  on  the  strain  hardening 
exponent  and  mildly  dependent  on  the  ratio  aQ/E.  On  the  basis 
of  Equations  6  and  8,  <5^,  the  crack  tip  opening  displacement,  is 
also  a  parameter  which  characterizes  the  intensity  of  the  stress- 
strain  HRR  field.  The  above  HRR  field  equations  are  applicable 
only  for  the  case  of  stationary  cracks. 

Much  of  the  current  research  in  modeling  "High-FCGR" 
behavior  utilizes  the  assumption  that  the  HRR  field  equations  are 
applicable  to  quasi-static  moving  cracks.  Paris13  has  suggested 
an  argument  to  justify  the  J-Integral  for  FCGR  studies  based  on 
the  following  logic: 

For  high  rates  of  growth,  the  crack  tip  moves 
ahead  during  each  cycle  into  relatively  new  material  in 
terms  of  plastic  deformation,  compared  to  the  intense 
deformation  it  will  sustain  at  the  crack  tip  during  the 
next  cycle.  Thus,  during  the  next  cycle,  past  history 
will  not  be  significant  compared  to  the  loading,  which 
is  then  being  sustained.  As  long  as  a  moving  crack  is 
considered,  it  may  be  possible  to  neglect  past  history 
including  unloading  in  a  J-Integral  analysis  and 
characterization  of  material  behavior. 
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Rice  suggests  that  a  crack  tip  integral  similar 

to  J  can  be  defined  for  cyclic  loadings  if  the  alteration  Ao^ 

of  stress  (i.e.  the  stress  range)  following  any  load  maximum  or 

minimum  is  related  to  the  corresponding  alteration  of  strain  Ae^, 

for  those  alterations  actually  experienced  by  the  material,  in 

such  a  way  that  /Ao^  dAe^j  is  a  function  only  of  the  Ac's 

and  not  of  position  in  the  material.  To  date,  no  publication  has 

appeared  in  which  the  AJ-Integral  has  been  calculated  by  applying 

the  integral  formulation  presented  in  Equation  7  to  the  crack 

tip  region  for  cyclic  problems.  Normally,  the  approach  taken  for 

such  loadings  is  to  use  the  potential  energy  property  of  the 

J-Integral  to  calculate  what  is  termed  an  operational  definition 

of  the  J-Integral.  In  a  review  of  Dowling  and  Begleys'  results 
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obtained  on  compact  type  specimens,  Parks  observed  that  objections 
to  the  operational  use  of  the  parameter  AJ  can  be  rationalized 
away  if  one  notes  the  fact  that  stable  hysteresis  loops  (macro¬ 
scopic  load/displacements)  are  nearly  symmetrical  for  reversed 
stressing.  This  observation  implies  that  most  material  points 
are  seeing  fully  reversed  (zero  mean)  deviator ic  stressing. 

The  stable  hysteresis  loop  hypothesis  leads  to  a  strain  energy 

3W* 

fvxnction  (W*)  approach  such  that  Aa  *  y  (Acj  ~  Ac  (Ae_)  exists  for 
all  material  points  in  the  crack  tip  region. 

2.2  CRACK  TIP  APPROACH 

As  discussed  earlier,  there  exists  an  intensely  deformed 
region  at  the  crack  tip.  The  deformation  processes  in  this 
intensely  deformed  process  zone  are  believed  to  control  the  crack 
growth  behavior.  Instead  of  correlating  the  crack  growth  with 
the  parameters  that  define  the  stress  field  surrounding  the 
process  zone,  it  should  be  possible  to  correlate  the  crack  growth 
with  parameters  identified  within  the  process  zone  directly. 

This  section  identifies  currently  available  crack  tip 
parameters.  Table  2  presents  some  of  these  parameters  which  have 
been  described  in  the  literature. 
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TABLE  2 

FRACTURE  MECHANICS  CRACK  TIP  PARAMETERS 


Crack  Tip  Parameters 
Near  Tip  Strain 

Crack  Tip  Opening  Displacement  (CTOD) 
Crack  Tip  Opening  Angle  (CTOA) 

Crack  Opening  Angle  (COA) 

Crack  Tip  Nodal  Force 

Energy-Based  Parameters 

Crack  Separation  Energy  Rate 

Change  in  Energy  Rate  in  Process  Zone 

Generalized  Energy  Release  Rate  at  the 
Process  Zone  (G  =  G,  +  G  ) 

A  p 


2.2.1  Crack  Tip  Parameters 

A  series  of  investigations  have  attempted  to  develop 
a  cause-effect  model  for  crack  advance  by  correlating  the  level  of 
stress-strain  displacement  related  parameters  in  the  crack  tip  region 
with  the  amount  of  crack  extension.  In  this  subparagraph,  those 
investigations  most  relevant  to  the  fatigue  crack  growth  problem 
are  considered. 

2 . 2 . 1 . 1  Near  Tip  Strain  Parameter 

To  enhance  our  understanding  of  fatigue  crack 

growth  behavior,  several  investigators  have  employed  a  critical 

strain  model.  The  critical  strain  model  is  based  on  the  concept 

that  an  element  of  material,  which  is  one  da/dN  wide,  will  fail  when 

the  strain  associated  with  that  location  in  the  material  reaches 

a  critical  level.  The  various  critical  strain  models  are  based  on 
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the  level  of  accumulated  plastic  strain  ,  maximum  strain  '  , 

or  strain  ranges  associated  with  inducing  a  sufficient  amount  of 
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fatigue  damage  m  the  element  .  The  strain  range  modeling 
approach  might  be  more  properly  called  a  fatigue  damage  approach 
because  failure  in  the  material  element  is  hypothesized  after  a 
sufficient  amount  of  fatigue  damage  has  been  induced  by  the  intense 
cyclic  action  as  the  crack  propagates  toward  the  material  element. 

The  results  of  this  approach  have  provided  better  correlations  with 
FCGR  damage  as  the  analytic  descriptions  of  the  crack  tip  stress 
and  strain  fields  have  improved. 

One  important  and  particularly  interesting 
feature  of  the  near  tip  strain  parameter  approach  is  that  the 
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strain  in  the  crack  tip  region  can  be  experimentally  measured  '  ' 

making  direct  correlations  between  this  parameter  and  the  resulting 

29 

FCGR  damage  possible.  Further,  as  Newman  has  shown,  a  maximum 
strain  criterion  can  be  used  as  a  simple  criterion  for  crack 
advance  in  a  cyclic  finite  element  program.  Newman  suggests  that 
when  the  maximum  strain  in  the  crack  tip  element  reaches  a  critical 
level  the  crack  tip  node  is  released  and  the  crack  advances  to 
the  next  node  in  line.  Using  this  approach,  Newman  has  been  able 
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to  successfully  reproduce  the  experimentally  observed  crack 

35 

opening  and  closing  stress  results  of  Elber 

2 . 2 . 1 . 2  Crack  Tip  Opening  Displacements  and  Angles 

The  near-tip  field  in  crack-growth 

situations  is  different  than  in  the  stationary  case.  There  is  no 

complete  description  of  the  stress  and  strain  fields  ahead  of  a 

36 

extending  crack.  Some  insight  has  emerged  due  to  work  by  Rice  , 
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Chitaley  and  McClintock  and  Amazigo  and  Hutchinson  .  These 

studies  revealed  a  milder  strain  singularity  of  the  form  £n  (1/r) 

3  6 

for  elastic-plastic  materials.  The  studies  by  Rice  based  on  J2 
flow  theory  of  plasticity  for  an  ideally  plastic  material  showed 
that  the  incremental  strains  at  the  crack  tip  point  are  related 
to  an  increase  of  the  crack  opening  displacement  d6,  and  the 
increment  of  the  crack  extension  da.  Hence,  the  strains  at  the 
crack  tip  are  uniquely  characterized  by  the  crack  opening  angle, 
d6/da. 

The  crack  tip  opening  displacement  ( <5 1 ) 
and  crack  tip  opening  angle  (at)  parameters  have  been  used  by 
Shih  et  al.^  and  Kanninen  et  al.^  to  predict  the  initiation  and 
stable  growth  of  cracks  subjected  to  monotonic  loading.  The  real 
interest  in  these  parameters  is  created  by  the  potential  for 
quantifying  the  level  of  crack  opening  either  from  experimental 
measurements  or  by  numerical  analysis.  The  concept  of  using  6^ 
and  at  as  FCGR  correlation  parameters  is  attractive  since,  for  each 
subcritical  crack  advance  (da/dN)  that  occurs  during  the  application 
of  a  cyclic  load,  the  crack  will  only  open  a  finite  amount.  One 
of  the  major  difficulties  that  restricts  the  utility  of  these 
parameters  is  the  lack  of  agreement  on  where  the  measurements 
should  be  made. 

Figure  12  provides  a  schematic  which 
describes  the  locations  where  measurements  should  be  taken  for 
comparison  to  available  analytical  results.  Two  definitions  of 
the  crack  opening  displacement  and  crack  opening  angle  have  been 
used  in  the  literature.  Crack  opening  displacement  (6^)  can  be 
measured  either  at  the  original  crack  tip  (6Q)  or  at  a  fixed 
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INITIAL 
CRACK  TIP 


FIXED  DISTANCE  FROM 
THE  CRACK  TIP 


Figure  12.  Crack  Tip  Displacement  Related  Parameters. 

Two  Definitions  of  Crack  Opening  Angle  and 
Displacement. 


angle  (45°)  from  the  crack  tip  (6^) .  The  two  possibilities  of  the 
crack  opening  angle  are  that  it  is  measured  either  at  the  original 
crack  tip  (a  )  or  at  a  fixed  distance  away  from  the  crack  tip  (ou. 

For  stable  crack  growth,  the  crack  opening 

39 

angle  stays  constant  as  is  observed  by  Shih  et  al.  and 
40 

Kanninen  et  al.  in  their  numerical  simulation  and  experimental 
measurements  of  a.  Therefore,  there  exists  a  theoretical,  numerical 
and  experimental  basis  for  the  use  of  a  crack  opening  angle  as  a 
crack  growth  criterion.  Also,  the  crack  tip  opening  angle  can 
be  measured  experimentally,  though  such  measurements  may  be 
difficult  at  high  temperatures. 

One  of  the  disadvantages  of  this  method  is 

the  lack  of  a  unique  definition  for  crack  opening  angle  or 

displacement.  Crack-tip  opening  angles  and  or  displacements  are 

defined  in  the  literature  at  a  fixed  distance  behind  the  crack-tip 

or  at  the  original  crack  tip.  The  theoretical  basis  for  using 

crack  tip  displacement  parameters  which  characterize  the  crack 

growth  is  derived  only  for  the  case  of  monotonically  increasing 

load,  but  this  solution  is  shown  to  be  valid  experimentally  and 

numerically  for  small  elastic  unloading  occuring  due  to  crack  growth 

in  the  process  zone.  Although  the  requirement  exists  for  more 

sophisticated  solutions  for  a  growing  crack  under  fatigue  conditions, 
39  4  0 

investigators  '  suggest  that  the  crack  opening  angle  may  correlate 
the  FCGR  behavior  generated  under  elastic-plastic  conditions. 

Finite  element  computational  procedures  have 
shown  that  the  COA  is  not  appreciably  sensitive  to  the  mesh  size. 
Further  work  is  necessary  in  this  area  to  see  whether  this  parameter 
is  applicable  for  FCG  situations. 
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2 . 2 . 1 . 3  Crack-Tip  Node  Force  (F) 


40 

Kanninen  et  al.  have  proposed  that  at  a 
constant  critical  force  at  the  crack-tip  node,  crack  growth 
occurs.  They  have  shown  that  this  critical  node  force  criterion 
performs  well  in  the  case  of  stable  crack  growth  situations  under 
elastic-plastic  conditions. 

The  main  disadvantage  of  this  method  is 
that  this  critical  force  is  sensitive  to  the  finite-element  mesh 
and  the  increment  step,  although  it  stays  constant  during  stable 
crack  growth  for  a  particular  mesh. 

Crack  tip  node  force  (F)  can  be  evaluated 
from  a  numerical  simulation  of  the  crack  growth.  Hence,  this 
parameter  could  be  extended  for  fatigue  crack  growth  situations. 

2.2.2  Energy-Based  Parameters 

The  generalization  of  Griffith's  energy  balance 

concept  to  account  for  nonlinear  material  behavior  has  been 

addressed  by  many  investigators  in  an  attempt  to  define  crack 

growth  parameters  (especially  stable  crack  growth  parameters) . 

An  early  continuum  approach  to  generalizing  the  Griffith  energy 
.  41 

balance  is  given  by  Rice  under  the  assumption  of  small  geometry 
changes.  The  resulting  energy  balance  takes  the  form  of 

G-P  =  2y  (  9  ) 

where  G  is  the  energy  release  rate,  P  is  the  plastic  work  rate, 
and  the  surface  energy  y  is  interpreted  as  the  work  per  unit  area 

required  to  create  new  surfaces.  For  materials  in  which  the 

.  .  .  41 

stress  saturates  to  a  finite  value  at  large  strain,  Rice  showed 

that  there  is  no  energy  surplus  from  the  continuum  calculation 

which  can  be  equated  to  the  work  of  separation.  At  the  other 

extreme,  when  the  material  is  elastic,  the  plastic  work  rate  term 

is  zero.  Therefore,  the  difference  G-P  varies  from  G  to  zero,  when 

the  hardening  behavior  of  material  varies  from  elastic  to  perfect 

plastic. 
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Kfouri  and  Miller  '  have  presented  finite 

element  solutions  for  the  crack  separation  energy  rate  (G^)  which 

is  related  to  the  work  done  in  separating  the  crack.  The  crack 

separation  energy  rate  (G^)  is  obtained  by  quasi-static  unloading  of 

a  finite  element  of  width  Aa  at  the  crack  tip,  while  the  applied 

stress  is  maintained  constant  on  the  structure.  Through  analyses 

of  the  relations  between  G^  and  Aa,  a  fracture  criterion  can  be 

suggested  by  taking  Aa  and  G^  as  characteristic  of  the  material. 

The  results  indicate  that  a  criterion  of  this  kind  predicts  stable 

crack  growth  when  the  step  size  is  small  compared  to  the  yield 

zone,  at  least  under  small  scale  yielding  conditions.  The  results 

also  show  that  unstable  fracture  occurs  immediately  when  the  step 

4  2  4  3 

size  is  greater  than  the  maximum  radius  of  the  plastic  zone  ' 
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Kanninen  et  al .  have  also  proposed  that  a  process 
zone  energy  release  rate  (G  ) ,  i.e.  the  change  in  energy  contained 
in  the  process  zone,  will  correlate  with  stable  crack  growth 
behavior  in  elastic-plastic  fracture.  The  crack  tip  process  zone 
is  the  region  at  the  crack  tip  where  non-proportional  loading 
occurs  during  crack  growth.  Some  believe  that  the  energy  flow  to 
this  process  zone  is  a  measure  of  the  strains  and  stresses  at  the 
extreme  tip  of  the  crack  and  hence  can  be  used  to  correlate  the 
fracture  process  occurring  at  the  crack  tip. 

The  size  of  the  process  zone  used  for  computing  the 

energy  change  of  this  zone  due  to  the  growth  of  the  crack  is  dependent 

on  the  interpretation  of  the  investigators  who  have  studied  this 
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approach.  Kanninen  et  al.  define  the  size  of  the  process  zone  as 

equal  to  the  plate  thickness  for  their  calculations,  whereas  Shih 
3  9 

et  al.  '  define  its  size  as  five  times  the  crack  tip  opening  displace- 
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ment  (56^) .  Green  and  Knott  define  a  process  zone  by  considering 
the  spacing  of  inclusions  and  ductile  fracture  by  cavities.  In 
their  process  zone,  strains  are  large  and  approximately  equal  to 
the  fracture  strain.  For  materials  with  a  low  work  hardening 
rate,  they  find  that  the  size  of  the  process  zone  is  twice  the 
crack-tip  opening  size;  for  a  high  work  hardening  rate,  the 
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material's  process  zone  size  is  less  than  twice  the  crack  tip 
opening  displacement.  Hence,  the  size  of  the  process  zone 
appears  to  be  a  material  parameter  which  depends  on  inclusion 
spacing. 

During  stable  crack  growth,  remains  constant, 
so  G  may  be  a  reasonable  parameter  to  characterize  the  crack 

tr 

growth.  Energy  criteria  are  evaluated  numerically  and  the  major 
disadvantage  of  this  parameter  is  that  it  is  sensitive  to  the 
mesh  size. 


as 
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A  generalized  energy  release  rate  (G)  is  defined 


G  = 


(10) 


which  is  equal  to  the  energy  flow  across  the  boundaries  of  the 
process  zone  during  the  elastic-plastic  stable  crack  growth  (See 
Figure  13) .  For  stable  crack  growth,  this  value  of  G  seems  to 
remain  constant.  The  computed  value  of  G  is  insensitive  to  the 
finite  element  mesh  size  and  thus  this  parameter  may  provide  a 
reasonable  criterion  for  studying  elastic-plastic  fatigue  crack 
growth  behavior. 


2.3  EVALUATION  OF  EPFM  PARAMETERS 


In  order  to  correlate  high  fatigue  crack  growth  rates  using 
EPFM  parameters,  it.  is  necessary  to  have  the  capability  to  calculate 
these  parameters.  The  results  of  continuum  analyses  of  elastic- 
plastic  crack  tip  material  behavior  are  limited  to  idealized 
geometries,  material  behavior,  and  loading  conditions.  Even  in 
those  limiting  cases,  a  complete  quantification  of  the  behavior 
by  analytical  methods  is  practically  impossible.  Among  a  number 
of  approximate  numerical  analysis  techniques,  the  finite  element 
analysis  provides  the  best  approach  to  this  task.  In  this  para¬ 
graph,  some  important  considerations  associated  with  applying  the 
finite  element  method  for  crack  tip  analysis  are  reviewed.  The 
results  presented  in  this  section  deal  with  monotonic  loading 
conditions  since  cyclic  crack  tip  analysis  is  still  in  a  state  of 
infancy.  We  expect  that  cyclic  loading  will  further  complicate 
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Figure  13.  Generalized  Energy  Release  Rate  (G) 


finite  element  analysis  of  the  crack  tip  beyond  that  which  is 
illustrated  for  monotonic  loading. 

2.3.1  Crack  Tip  Analysis  by  Finite  Elem ^nts 

In  the  modeling  of  an  elastic-plasLic  continuum  which 
contains  a  crack,  two  features  of  the  crack  geometry  must  be  given 
special  consideration:  (1)  the  crack  tip  singularlity ,  and 
(2)  the  extent  of  crack  tip  shape  changes  (blunting) .  In  the 
paragraphs  below,  we  will  outline  several  finite  element  modeling 
approaches  which  deal  with  these  geometric  features. 

2.3. 1.1  Crack  Tip  Singularity  Modeling 

For  problems  with  singularities,  Tong  and 
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Pian  have  shown  that  the  convergence  rates  of  the  finite  element 

solutions  are  dominated  by  the  nature  of  singularity  near  the  crack 

tip.  Unless  these  singularities  are  properly  represented  in  the 

assumed  solution  functions,  even  with  regular  so-called  high-accuracy 

elements  (with  higher  order  polynomials  as  interpolates) ,  the  rate 

of  convergence  of  solutions  cannot  be  improved.  However,  when  the 

functions  for  the  singularity  are  incorporated  explicitly  in  the 

assumed  solution  in  the  near  crack  tip  elements,  in  general  the 

satisfaction  of  the  convergence  criteria  for  the  finite  element 
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solution  is  not  guaranteed.  Current  finite  element  approaches 
for  solving  crack  tip  problems  fall  into  two  general  classes. 

The  first  approach  ignores  the  singular 

nature  of  the  solution  (Equation  5)  in  those  finite  elements  near  the 
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crack  tip  and  uses  regular  assumed  solution  functions  .  This 

approach  ensures  that  the  finite  element  convergence  criteria  are 
satisfied  but  requires  a  very  fine  finite  element  grid  near  the 
crack  tip  to  obtain  reasonable  results  at  or  close  to  the  crack  tip. 

A  disadvantage  of  this  approach  is  that  it  is  costly. 

The  second  approach  incorporates  appropriate 
.  ci.cq 

functions  for  the  singularity  in  the  crack  tip  elements 

in  one  of  the  following  ways: 
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(a)  Exact  crack  tip  singularities  are 
included  in  the  near  crack  tip  elements.  In  this  way,  the 
singular  solution  can  be  represented  accurately  witn  relatively 
few  elements,  but  the  compatibility  of  displacements  between 
crack  tip  and  regular  elements  is  not  assured. 

(b)  Crack  tip  elements  are  represented 
by  degenerated  isoparametric  elements  or  quarter  point  8-noded 
isoparametric  quardilateral  elements.  This  approach  requires  no 
special  coding  for  elements  at  the  crack  tip,  introduces  no 
incompatibility,  and  involves  no  extra  cost. 

(c)  The  third  approach  utilizes  a  hybrid 
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finite  element  in  the  modeling.  The  solution  functions  in 

these  crack  tip  elements  are  chosen  so  that  the  singularity  is 
appropriately  modeled.  The  interelemcnt  continuity  condition, 
which  is  required  as  a  convergence  criterion,  is  not  violated 
due  to  the  presence  of  the  singular  functions.  The  interelement 
continuity  condition  is  enforced  as  an  a  posteriori  condition  of 
extremization  of  the  modified  energy  function  by  a  Langrange 
multiplier  method. 

Table  3  summarizes  the  advantages  and 
disadvantages  of  the  two  general  classes  of  finite  element  approaches. 

2. 3. 1.2  Crack  Tip  Blunting 

The  other  problem  involved  in  the  finite 
element  formulation  for  elastic-plastic  crack  analyses  is  the 
development  of  a  finite  element  grid  which  will  model  the  configura¬ 
tion  of  the  crack  tip  blunting  during  loading.  Large  deformation 
analysis  is  required  in  the  investigation  of  the  crack  tip  blunting. 
McMeeking  studied  the  crack  tip  finite  deformation  configuration 
for  small  scale  yielding  by  using  the  very  fine  mesh  of  4-node 
quadrilateral  isoparametric  finite  elements  shown  in  Figure  14. 

After  a  sufficient  amount  of  load  was  applied,  he  expected  that 
tne  steady-state  solution  for  contained  yielding  to  be  largely 
independent  of  original  geometry  of  the  crack  tip.  When  the  crack 
tip  opening  was  more  than  double  the  original  opening,  the 
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COMPARISON  BETWEEN  FINITE  ELEMENT  MODELS  FOR  FRACTURE  ANALYSIS 


Figure  14.  Undeformed  Configuration  of  Finite  Element 
Mesh  for  the  Notch-Blunting  Solution 
(Reference  58). 

solutions  for  continuing  load  increments  settled  down  to  a  steady 
level  in  which  increments  of  all  quantities  with  length  dimensions 
(such  as  crack  tip  opening  and  plastic  zone  size)  are  scaled  with 
increments  of  J/a0.  McMeeking  computed  the  J  parameter  on 
contours  remote  from  the  tip.  Thus,  above  a  certain  load  level, 
McMeeking' s  finite  element  grid  appears  to  represent,  with 
reasonable  accuracy,  the  blunting  behavior  of  the  initially  sharp 
crack  tip. 
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Shih  et  al.  utilized  degenerate  8-node 
isoparametric  quadrilateral  elements  as  discussed  by  Barsoum^ 
for  the  singular  crack  tip  element  as  shown  in  Figure  15.  In 
the  figure,  only  the  corner  nodes  are  indicated.  The  crack  tip 
elements  with  triangular  shape  in  Figure  15a  are  8-node  quadri¬ 
lateral  elements  in  which  the  two  corner  nodes  and  the  corresponding 
mid-side  node  are  initially  collapsed  to  a  common  point  at  the 
crack  tip.  As  the  load  is  increased,  the  crack  tip  blunting  is 
modeled  by  separating  the  nodes  at  the  common  points,  as  shown  in 
Figure  15b. 


Figure  15.  Finite  Element  Mesh  for  Crack-Tip  Blunting 
(Reference  39)  . 


2.3.2  Numerical  Estimation  of  the  J-Integral 

Shih  et  al.22  have  proposed  a  procedure  to  estimate 
crack  tip  parameters  J  and  which  were  generated  for  several 
cracked  structures  manufactured  from  materials  that  can  be 
modeled  using 
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The  structures  considered  were  center-cracked,  single-edge-notched, 
compact  tension  and  three  point  bend  specimen  geometries.  The 
Shih,et  al.  procedure  is  based  on  a  series  of  solutions  to 
boundary  valued  problems  where  only  the  strain  hardening  part, 
a  (a/a0)n  ,  of  Equation  11  is  used  to  model  the  material.  They 
have  observed  that  when  the  solutions  involve  a  monotonically 
increasing  single  load  parameter,  P,  all  stress  and  displacement 
parameters  increase  in  direct  proportion  to  the  load  parameter 
raised  to  a  power  dependent  on  the  strain  hardening  exponent,  n. 

The  value  of  J  obtained  for  a  center-cracked  panel  made  from  a 
strain  hardening  material  is  given  by: 
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J  =  a  °o  eo  a  (  TT  >  '  9i<5'n>*<^)n  1  (12) 

where  a  is  the  half  crack  length,  b  is  the  half  width  of  the 
specimen,  is  a  function  of  ^  and  n.  Equation  (12)  will  provide 
a  simple  functional  relationship  which  can  be  used  to  compute  J 
for  different  load  levels.  These  simple  functional  relationships 
of  J  or  for  different  geometrical  configurations  would  provide 
a  simple  method  of  evaluating  these  parameters  for  different  crack 
lengths  and  materials  (different  values  of  a  and  n) . 

Once  the  solutions  given  by  Equation  12  are 
established,  these  can  be  employed  together  with  linear  elastic 
solutions  to  produce  approximate  formulas  for  the  estimation  of 
J  and  crack  opening  displacement  which  are  valid  in  the  transition 
region  between  small-scale  yielding  and  large-scale  yielding. 

2-3.2  Concluding  Remarks  on  Finite  Element  Results 

Numerous  reports  on  the  finite  element  analyses  of 
elastic-plastic  fracture  problems  are  presently  available.  The 
potential  of  the  finite  element  method  as  a  tool  to  obtain 
quantitative  solutions  of  practical  problems  with  complicated 
geometries,  and  to  verify  approximate  continuum  analysis  results 
numerically,  is  considerable. 

The  repeatability  of  the  finite  element  solutions 
i-or  the  linear  elastic  fracture  analysis  using  different  analytical 
techniques  is  quite  good.  However,  for  elastic-plastic  fracture 
analysis,  different  finite  element  methods  do  not  always  provide 
accurate  and  unambiguous  solutions. 

This  fact  is  evident  from  the  work  of  Wilson  and 
Osias^  which  describes  the  results  of  a  ASTM  round-robin 

vest  tion  which  was  conducted  (around  1973-1974)  to  compare 
the  results  of  ten  different  state-of-the-art  finite  element 
solutions  of  a  cracked  three-point  bend  specimen.  As  can  be 
seen  from  Figure  16,  the  discrepancies  between  the  various  analyses 
at  high*  Load  levels  are  significant. 
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The  letters  on  each  curve  of  Figure  16  represent 
the  analyses  described  in  Table  4.  Unfortunately,  no  experimental 
results  are  available  for  comparison.  ASTM  Subcommittee  E24.08 
on  Elastic  Plastic  Fracture  is  currently  involved  in  a  new 
round-robin  investigation  of  Finite  Element  Computational 
Capability. 
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Figure  16. 


Crack  Mouth  Opening  Displacements  as  a 
Function  of  Magnitude  of  Applied  Load 
(Reference  60) . 
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SECTION  III 

SUMMARY  OF  AVAILABLE  DATA 

In  this  section  of  the  report,  FCGR  data  are  identified 
that  have  been  collected  under  conditions  which  induce  extensive 
crack  tip  yielding.  The  section  is  organized  so  that  attempts  to 
correlate  FCGR  data  with  LEFM  parameters  are  presented  first. 

These  correlations  are  followed  with  a  summary  of  currently  available 
J-Integral  work.  Finally,  all  other  known  correlations  of  FCGR  data 
with  other  parameters  are  summarized. 

3.1  THE  FCGR  LEFM  PARAMETER  CORRELATIONS 

3.1.1  Stress  Intensity  Factor  Correlations 

4  5 

As  indicated  in  Section  1,  Frost  et  al.  ,  Dubensky  , 
and  Dowling  and  Walker^  developed  FCGR  correlations  based  on  the 
LEFM  parameter  AK.  Each  investigator  found  that  some  ductile 
materials  in  the  high  FCGR  regime  do  not  exhibit  a  correlation 
that  is  independent  of  the  maximum  applied  stress.  Dubensky 
attempted  to  correct  for  the  effect  of  plasticity  by  modifying 
the  stress  intensity  factor  calculation  through  the  use  of  an 
effective  crack  length.  Dubensky' s  effective  crack  length  was 
obtained  by  extending  the  observed  crack  length  values  by  a 
calculated  plastic  zone  size. 

The  2024-T3  aluminum  data  which  exhibited  the 
maximum  stress  influence  identified  in  Figure  4 ,  still  showed 
the  stress  level  influence  even  with  Irwin,  Dugdale  and  Newman 
type  plasticity  modifications  considered.  Dubenski's  data  provide 
a  valuable  base  for  calculations  and  the  details  of  his  test 
program  are  summarized  in  Table  5. 


In  Dubensky's  test  program,  FCGR  data  were  collected 
in  12  inch  wide  center-cracked  panels  of  0.090  inch  thickness 
for  high  net  section  stress  conditions.  Initial  starting  half- 
crack  length  sizes  were  in  the  vicinity  of  0.080  inch  for  all 
tests.  Described  in  the  Table  5  is  a  summary  of  the  initial 
and  failure  net  section  stress  conditions  for  the  FCG  tests 
conducted  on  the  two  materials.  Also  provided  are  fracture  toughness 
calculations. 

TABLE  5* 

SUMMARY  OF  DUBENSKY'S  TEST  RESULTS 


Material 

Stress  Ratios 
Considered 

°nat  | 

oy#  1 a«0.08  inch 

a 

net. 

5  1 a-a - 

ys  c 

(Average) 

Average 

( inch) 

Average 

*C 

Ksi/in. 

7076-T6 

0 

0.74 

0.77 

0.345 

57.3 

0,0.33,0.3,0.7 

0.90 

0.833 

0.258 

53.8 

0,0.33,0.5,0.7 

0.87 

0.385 

0.199 

51.3 

0,0.33,0.5,0.7 

0.94 

0.953 

0.149 

47.8 

0,0.33,0.5 

0.97 

0.976 

0.113 

45.5 

2024-T3 

0,0.33,0.5,0.7 

0.59 

1.18 

3.02 

110. 

0,0.33,0.5,0.7 

0.79 

1.08 

1.67 

96.2 

0,0.33,0.5,0.7 

0.99 

1.06 

0.460 

60.6 

0,0.33,0.5,0.7 

1.04 

1.06 

0.196 

41.1 

*  1  in  =  25.4  mm;  lKsi/Tn  =  1.1  MN/m^ 


As  can  be  noted,  net  section  stresses  at  failure  are  above  75 
percent  of  the  yield  strength  for  the  7075  alloy  and  above  105 

for  the  2024  alloy.  Fracture  toughness  is  noted  to  vary 
with  the  fracture  half-crack  length  (a^)  for  both  materials  but 
more  dramatically  for  the  2024  alloy.  The  crack  length  versus 
cycles  data, on  which  the  FCGR  data  were  generated  ,  are  contained 
in  the  report.  The  2024-T3  aluminum  alloy  data  for  the  three 
highest  maximum  stress  conditions  are  presented  in  Figure  XI  as 
a  function  of  As  can  be  seen  from  the  Figure,  the  normal 

stress  ratio  (R)  effect  disappears  as  the  maximum  stress  and 
crack  growth  rate  increase.  For  the  two  highest  maximum  stress 
conditions,  the  FCGR  cluster  about  a  single  trend  line  which  is  a 
characteristic  for  that  maximum  stress  condition.  Figure  17 
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Highest 


provides  direct  evidence  that  it  may  be  necessary  to  describe 
FCGR  above  a  given  level  using  a  two  parameter  approach  where  one 
parameter  is  a  fracture  and  the  other  is  a  strength-related 
parameter. 

Wilhem  and  Ratwani  have  been  actively  conducting 

fracture  resistance  testing  for  the  last  several  years**"*" '  . 

Recently  they  attempted**^  to  correlate  the  FCGR  data  collected 
5 

by  Dubensky  on  2024-T3  aluminum  with  some  of  their  monotonic¬ 
loading,  fracture-resistance-curve  test  results.  Figure  18  shows 
that  they  have  had  reasonable  success  in  correlating  the  tearing 
crack  movement  (Aa)  under  monotonic  loading  conditions  with  the 
average  crack  advance  per  cycle  (Aa=da/dN)  for  fatigue  loading 
conditions.  For  both  types  of  loadings,  the  maximum  stress 
intensity  factor  (K  =K  )  was  used  to  portray  the  results. 

While  Wilhem  and  Ratwani ' s  results  appear  to  be  at  variance  with 
the  data  presented  in  Figure  17,  it  must  be  pointed  out  that 
Wilhem  and  Ratwani  used  Dubensky' s  lower  maximum  (gross)  stress 
level  FCGR  data  which  appear  to  correlate  reasonable  well  using 
LEFM  correlations  (See  Figure  4) . 

Because  of  somewhat  striking  differences  between 
the  FCGR  behavior  generated  under  low  and  high  gross  stress 
conditions,  one  might  first  attempt  a  correlation  on  the  basis  of 
net  section  stress  levels.  As  can  be  noted  from  the  Table 
summarizing  Dubensky' s  data,  such  an  approach  does  not  appear  too 
promising.  For  purposes  of  establishing  some  criterion  that  will 
describe  the  boundary  between  when  FCGR-LEFM  correlations  can  or 
can  not  be  used  Frost  et  al.4  and  Dowling®6  suggested  zone  criteria. 

4 

The  Frost  et  al.  criterion  (from  Equation  2) : 


r 


Y 


1  ,  Kmax  s  a 

Tff  1  -  7 


(13) 


for  the  Ixmiting  plastic  zone  size  is  suggested  as  an  upper  bound 
for  restricting  the  utility  of  the  LEFM  parameters.  For  the  most 
part,  Dubensky' s  high  FCGR  data  generated  for  the  lower  (gross) 
stress  maximum  levels  in  the  2024-T3  aluminum  exhibit  higher  (r^/a) 
levels  than  suggested  by  Equation  13;  conversely,  the  FCGR 
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data  generated  for  the  higher  maximum  stress  levels  exhibit  levels 
of  the  (r^/a)  parameter  more  nearly  in  line  with  the  Frost  et  al.^ 
criterion. 

66 

Dowling  more  recently  considered  the  problem  where 
geometric  dimensions  as  a  whole  must  be  considered  in  estimating 
the  upper  limits  of  LEFM  applicability.  In  this  study,  compact 
specimens  were  used  and  the  specific  focus  was  on  the  ratio  of 
plastic  zone  size  to  remaining  ligament,  width  minus  crack  length 
(W-a) .  Dowling  showed  that  an  appropriate  criterion  which  restricted 
unbounded  plasticity  is  given  by 


1  ,  max  . 2  ^  W-a  , .  . . 

ry  ’  “ST  <  —  >  <  T  <14) 

rie  subsequently  employed  this  criterion  to  identify  FCGR  data, 
generated  for  Man  Ten  Steel,  which  could  not  be  correlated  utilizing 
the  AK  (LEFM)  parameter  6  (See  Figure  5) . 


3 . 1 . 2  Strain  Intensity  Factor 

6  7 

In  1963,  McEvily  and  Boettner  reported  that 
they  were  able  to  describe  FCGR  behavior  of  copper  between  0.002 
and  0.02  mm/cycle  (intermediate  range)  using  a  strain  intensity 
range  approach.  More  recently,  Kitagawa  et  al.oS  have  attempted 
to  correlate  the  FCGR  behavior  associated  with  small  micro-cracks 
with  both  a  plastic  strain  intensity  range  and  total  strain 
intensity  range  with  some  degree  of  success  using  a  mild  steel 
(SB-22)  in  the  FCGR  range  between  10  and  10  mm/cycle.  Their 
studies  are  basically  concentrated  on  low-range  of  fatigue  crack 
growth  rates,  and  hence  the  plastic  zone  size  may  have  been  small 
compared  to  the  surrounding  elastic  field. 

Strain  intensity  is  a  measure  of  the  intensification 
of  strains  at  the  crack  tip  (K£-  eB/ifa)  .  Th->  computational  schemes 
associated  with  LEFM  parameters  which  are  documented  in  handbooks 
are  used  to  find  the  geometrical  component  of  the  strain  intensity 
factor  where  the  strain  is  either  the  plastic  struin  range  or  the  total 
strain  range. 
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3.2 


J-INTEGRAL  PARAMETERS 


As  suggested  by  Equation  6 ,  the  magnitude  of  the  stresses 
and  strains  in  the  elastic-plastic,  near  crack  tip  stress  and 
strain  fields  can  be  related  to  the  value  of  the  J-Integral,  the 

strain  hardening  exponent  (N) ,  the  yield  strength  level  (o  )  and 

13  18  19  ^ 

the  yield  strain  level  (eQ)  '  '  .  The  J-Integraj.  is  a  path 

independent  line  integral  described  by  the  expression  presented 

20 

in  Equation  7.  In  1968,  Rice  published  a  paper  describing  the 
analytical  properties  of  the  parameter  wherein  it  could  be 
(a)  correlated  with  the  Energy  Release  Rate  (G)  for  linear  elastic 
structures,  and  (b)  utilized  with  a  nonlinear  elastic  material 
analysis  to  describe  the  response  exhibited  by  elastic-plastic 
materials.  These  special  properties  led  Begley  and  Landes 
to  propose  how  this  parameter  could  be  utilized  for  measuring  the 
fracture  toughness  of  ductile  materials.  In  particular,  Begley 
and  Landes  proposed  using  an  operational  definition  of  the 
J-integral  where  it  is  assumed  that  the  crack  does  not  extend  and 
that  there  is  no  material  unloading  within  the  path  of  the  integral. 

3.2.1  Operational  J--Techniques  -  Summary 

?5 

Following  Dowling  and  Begley’s  *  presentation 

at  a  fracture  mechanics  symposium  in  1S74,  this  parameter  has 

received  a  great  deal  of  attention.  Dowling  in  subsequent 

work  ^  has  been  able  to  show  that  his  parameter  (AJ) 

correlates  FCGR  behavior  developed  front  several  different  crack 

geometries  (including  microcrack  sizes)  that  have  been  subjected 

to  constant  stress  amplitude  cyclic  histories.  Figure  19 

describes  the  success  experienced  by  utilizing  this  particular 

parameter.  The  reader  will  note  from  the  figure  that  in  the 

overlapping  range,  one  can  not  distinguish  differences  between 

those  data  derived  from  "AJ"  tests,  where  extensive  plastic 

deformation  is  observed,  and  from  tests  in  which  the  plasticity 

is  contained  under  small  scale  yielding  conditions.  Others  have 

duplicated  the  general  observations  of  Dowling.  In  particular, 

74 

Sadananda  and  Shahinian  have  recently  shown  its  application  to 
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<Sa/dN.  Crack  Growth  Rate,  inches/ cycle 


AJor  l AK> Z/ E.  in-lb/ in2 


Figure  19.  Fatigue  Crack  Growth  Rate  Versus  Cyclic 
J  for  Various  Geometries  (Reference  71) 


Udimet  700  at  850°C  for  a  number  of  different  constant  stress 
amplitude  histories. 

The  principal  advantage  of  the  AJ  approach  is 
that  it  appears  to  work  and  that  it  can  be  related  in  a  quasi¬ 
sense  to  the  plastic  stress-strain  field  at  the  crack  tip. 

A  major  disadvantage  of  the  approach  is  that  there 
are  several  computational  schemes  that  are  currently  available 
for  actually  calculating  numerical  values  for  this  parameter. 

75 

These  calculation  schemes  use,  as  their  basis,  the  Rivl in-Thomas 
compliance  technique  applicable  to  calculating  the  energy  release 
rate  (G)  for  nonlinear  elastic  material.  A  schematic  illustrating 
the  essence  of  the  compliance  technique  is  presented  in  Figure  20. 
The  difference  in  the  area  between  the  two  load-displacement 
curves  is  equal  to  the  energy  release  rate  for  nonlinear  elastic 
materials.  The  energy  release  rate  has  been  shown  by  Rice  20 
to  be  directly  equivalent  to  the  J-Integral.  While  the  calculation 
performed  under  monotonic  loading  is  rather  straightforward,  the 
analyst  must  define  a  common  origin  for  the  measured  cyclic  load- 
displacement  responses  before  starting  the  calculation.  Two  schemes 
have  been  suggested  for  defining  the  common  origin  for  the  cyclic 
(loading  portion)  load-displacement  curves  exhibited  by  compact 
specimens  under  constant  amplitude  loading;  these  are  presented 
in  Figure  21.  An  additional  disadvantage  that  might  also  be 
associated  with  this  seemingly  popular  technique  is  that  to  date 
all  calculations  of  AJ  have  been  based  on  experimentally  derived 
load-displacement  behavior . 

From  a  forward  looking  point  of  view,  the  experi¬ 
mental  measurement  of  AJ  represents  an  advantage,  since 
analytically  and  experimentally  derived  values  for  this  parameter 
can  be  compared. 
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Figure  20.  Procedure  for  Calculating  Operational  Values  of  the  J-Integral. 
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Figure  21.  A  Schematic  Description  of  the  Two  Approaches 
Used  to  Determine  Operational  Values  of  the 
J-Integral  -  After  Reference  74. 
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3.2.2  Operations 1-J  Data  Summary 


Dowling  and  Begley  pioneered  the  use  of  the 

Operational-J  approach  to  the  FCGR  correlation  in  the  presence  of 

gross  plasticity.  Compact  tension  fracture  specimens  of  A533B 

steel  were  subjected  to  gross  cyclic  plastic  deformations  and  FCGR 

up  to  0.01  in/cycle  were  obtained.  The  FCGR  data  correlates  with 

the  Operational  AJ  parameter  which  is  estimated  from  load  versus 

deflection  hysteresis  loops.  For  two  gross  plastic  specimens, 

AJ  vs.  N  and  a  vs.  N  plots  are  given  in  this  paper.  For  one  of 

the  above  specimens,  5  (maximum  deflection  of  the  COD  gauge) 

max 

vs.  N  cind  crack-opening  6  vs.  N  plots  are  also  given.  For  gross 
plastic  testing  conducted  under  deflection  control  to  a  sloping 
line,  FCGR  correlate  extremely  well  with  the  Operational  AJ  values. 
The  tests  conducted  with  load  control  to  a  sloping  line,  for  large 
AJ  values,  FCGR  deviate  drastically  from  the  straightline  behavior 
observed  in  da/dN  vs.  AJ  plots  derived  from  deflection  control 
tests.  Otherwise,  FCGR  for  both  K  tests  and  J  tests  is  on  the 
same  power  law  curve. 
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Dowling  continued  to  use  the  Operational-J 
approach  to  evaluate  the  effect  of  geometry  on  FCGR  behavior. 

Fatigue  crack  growth  data  were  obtained  for  center  cracked 
specimens  of  A533B  steel  subjected  to  elastic-plastic  cyclic  load¬ 
ing.  Cyclic  J-Integral  values  estimated  from  load  versus  deflection 
hysteresis  loops  are  correlated  with  these  growth  rate  data.  These 
results  are  in  agreement  with  elastic-plastic  data  on  CT  specimens 
and  also  with  linear  elastic  data  on  large  size  CT  specimens.  The 
tests  were  conducted  with  the  load-line  deflection  controlled  to  a 
sloping  line  on  load  deflection  planes.  Hysteresis  loops  and  crack 
lengths  are  also  given  for  two  elastic-plastic  specimens  by  this 
paper . 
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Extension  of  Operational-J  to  the  type  of  crack 

growth  that  occurred  during  low-cycle  fatigue  testing  of  smooth 

axially  loaded  specimens  (LCF)  fabricated  from  A533B  Steel  was 
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investigated  by  Dowling  .  Surface  crack  lengths  were  monitored 
periodically  using  cellulose  acetate  replicas.  The  estimation 
procedure  used  to  obtain  the  Operational-J  for  small  cracks  is 
summarized  in  Figure  22.  Experimental  crack  growth  for  different 
strain  levels  are  correlated  by  Operational-J  integral  as  shown 
in  Figure  23.  When  the  crack  length  is  small  (SL  <  0.007  in), 
the  resulting  FCGR  data  are  found  to  be  substantially  higher  than 
expected  on  the  basis  of  available  long  crack-low  stress  type  data 

When  the  estimation  procedure  for  calculating  AJ 

for  these  short  cracks  was  modified  so  that  crack  lengths  were 

8 

increased  by  an  effective  c rack  length  ( £  )  ,  El  Haddad  et  al. 
have  shown  that  the  crack  growth  data  from  the  small  cracks  in 
Figure  23  also  can  be  correlated  by  the  Operational-AJ  parameter 
as  shown  in  Figure  24.  The  value  of  £q  is  evaluated  from  the 
threshold  stress  intensity  and  the  endurance  limit  of  the  material 
El  Haddad  suggests  that  the  microstructural  features,  e.g.  £q 
value  is  a  material  constant  that  is  related  to  the  grain  size  for 
low  carbon  steel. 


Figure  22. 


Estimation  of  Small  Crack  AJ  from  Stress-Strain 
Hysteresis  Loops  (Reference  72) . 
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24.  Fatigue  Crack  Growth  Rates  as  a  Function 
of  Modified  AJ  (Reference  8) . 
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Figure  25.  Crack  Growth  Data  in  Terms 
Determined  Using  the  MerkL 
Approximation  (Reference  7 
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Sadananda  and  Shahiman  have  also  conducted  cyclic 

J-Integral  FCG  testing  on  Udimet  700  at  high  temperature  where 

the  possibility  of  a  larger  plastic  zone  is  likely.  Fatigue 

crack  growth  rates  were  analyzed  using  both  nonlinear  and  linear 

fracture  mechanics  where  the  tests  were  conducted  at  850°C 

under  a  load  control  mode.  Operational-AJ  was  estimated  both  by 

the  compliance  technique  illustrated  in  Figure  20  and  21  and  by 
7  6 

the  Rice  et  al.  estimation  .procedure  modified  by  the  Merkle  and 
77 

Corten  correction.  One  interesting  feature  of  the  Sadananda 
and  Shahinian  work  is  that  these  authors  present  their  FCGR  data 
in  terms  of  an  effective  stress  intensity  factor  (AK)  calculated 
using  the  expression: 

AK  (effective)  =  /AJ*E  (15) 

In  this  manner,  currently  available  FCGR  data  that  is  expressed 
in  an  LEFM  format  can  be  directly  compared  to  new  (higher  FCGR) 
data  without  converting  the  older  data  to  an  EPFM  format.  Shown 
in  Figure  25  are  FCGR  data  that  have  been  plotted  as  a  function 
of  the  effective  AK  parameter  calculated  using  the  Rice-Merkle- 
Corten  estimation  procedure.  The  figure  shows  that  these 
Udimet  700  data  are  described  reasonably  well  with  the  trend  line 
curve  established  using  the  more  tedious  compliance  technique 
(dotted  line  in  Figure) .  Furthermore,  when  the  AK  parameter 
is  calculated  using  LEFM  assumptions,  the  correlation  with  the 
EPFM  established  parameter-FCGR  data  values  is  poor. 

7  8 

These  authors  have  also  recently  shown  that  FCGR 
could  be  correlated  by  /AJ*E  for  cold  worked  Type  316  stainless 
steel  at  high  temperatures. 
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Mowbray  studied  fatigue  crack  growth  in  Cr-Mo-V  steel 

-4  -1 

in  the  high-rate  regime  of  2.5  x  10  to  2.5  x  10  mm/cycle 

with  a  compact  type  strip  specimen.  FCGR  data  were  correlated  as 

a  function  of  AJ  (or  AK2/E  for  elastic-plastic)  for  strip 

CT  specimens  both  in  elastic  and  elastic-plastic  tests.  Crack 

length  versus  number  of  cycle  curves  and  cyclic  hysteresis  loops 

are  presented  for  some  specimens.  For  elastic-plastic  high  growth 

rate  tests,  the  correlating  parameter  (AJ)  is  computed  using 

80 

the  estimation  procedure  suggested  by  Bucci  et  al. 
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As  a  consequence  of  specimen  geometry,  AJ  apparently 

seems  to  stay  constant  for  small  ranges  of  crack  growth.  However, 

a  careful  review  of  Mowbray's  work  reveals  that  the  specimen 

geometry,  material,  and  loading  combined  in  such  a  way  that 

(a)  the  LEFM  parameter  AK  was  essentially  constant  with  crack 

length,  and  (b)  the  nonlinear  plastic  deformation  was  limited. 

When  these  conclusions  are  reviewed  collectively,  the  values  of 

Mowbray's  Operational  AJ  are  within  twenty  percent  of  an  elastic 

2 

calculation  (AJ  =  (AK)  /E)  and  therefore  the  test  data  he  reports 
(as  constant  AJ)  are  really  the  result  of  a  constant  K  control 
condition. 
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Brose  and  Dowling  have  conducted  J-Integrai 
fatigue  testing  to  evaluate  the  effect  of  size  on  the  FCGR  of 
Type  304  Stainless  Steel.  The  da/dN  versus  AK  data  are  compared 
on  the  basis  of  several  size  criteria  which  are  intended  to  limit 
plasticity  and  hence  enable  linear-elastic  analysis  of  data. 

The  cyclic  J-Integral  method  of  testing  has  been  employed  in 
testing  specimens  undergoing  gross  plasticity.  The  data  from 
cyclic  J-Integral  tests  agree  well  with  the  linear-elastic  tests 
and  the  results  are  presented  using  da/dN  versus  AK  (/AJ*E)  plots. 
Hysteresis  loops  are  given  for  the  specimens,  which  were  tested 
using  the  cyclic  J-Integral  method  of  testing.  Load  line  deflection 
versus  crack  length  data  are  also  given  for  the  specimen  which 
was  tested  under  a  cyclic  J-Integral  method  of  testing. 
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Marschall  et  al.  recently  generated  fatigue  crack 
growth  for  a  cast  2-1/4  Cr-lMo  steel  at  an  elevated  temperature 
(540°C)  under  quasi-AJ  control  conditions.  A  typically  initial 
loading  cycle  is  shown  in  Figure  26.  Loading  was  reversed  at  C 
and  unloaded  to  the  point  E  such  that  the  maximum  compressive  load 
was  equal  to  the  prior  maximum  tensile  load.  The  point  D  is  an 
indication  of  the  beginning  of  the  crack  closure.  The  cyclic 
crack  growth  tests  were  conducted  with  a  maximum  compressive  load 
hold  which  corresponds  to  the  EF  in  the  load  displacement  diagram. 
The  specimen  is  loaded  from  F  to  H  such  that  the  AJ  evaluated 
from  the  area  under  GHN  is  equal  to  the  AJ  calculated  in  the 
previous  cycle.  This  process  has  been  continued  for  20  cycles 
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and  the  specimen  was  broken  open  to  show  the  crack  growth  process 
Marschall,  et  al.  results  indicate  that  the  crack  growth  rate 
results  from  the  tests  with  compressive  hold-time  greater  than 
one  day  can  be  correlated  by  the  Operational-AJ  parameter. 

They  also  reported  that  the  crack  growth  rate  may  not  be  constant 
over  the  entire  20  cycles  investigated,  but  may  diminish  as  the 
test  proceeds.  One  experimental  problem  experienced  by  these 
investigators  was  that  the  crack  only  advanced  in  the  interior, 
and  not  on  the  surface  of  their  compact  specimens. 


3.3  OTHER  FCGR  CORRELATIONS 


Only  two  investigations  have  been  identified  which  provide 

parameters  which  could  be  considered  as  alternates  to  tho 

Operational- J  Integral  parameter.  In  the  first  study,  high  FCGR 

behavior  for  a  mild  steel  was  described  as  a  function  of  the 
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crack  opening  displacement  .  The  data  resulting  from  this 
study  are  presented  in  Figure  27,  where  the  rate  of  fatigue  crack 
growth  is  roughly  proportional  to  (COD)  ** . 
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In  another  study,  Sadananda  and  Shahiman  attempted  to 
use  the  maximum  value  of  the  load-line  displacement  as  a 
correlation  parameter.  On  the  basis  of  the  FCGR  correlation 
described  in  Figure  28, it  appears  that  correlation  is  possible. 

These  authors  attempted  to  provide  a  correlation  with  this  parameter 
because  it  is  substantially  easier  to  measure  than  the  crack 
tip  opening  displacement,  especially  at  high  temperatures.  The 
authors  note  that  a  one-to-one  relationship  between  the  load- 
line  displacement  and  crack  tip  opening  displacement  may  not 
exist  when  the  plastic  zone  size  is  large.  Sadananda  and  Shahinian 
also  note  that  the  correlation  provided  by  the  load-line 
displacement  parameter  results  in  more  data  scatter  than  the 
correlations  resulting  from  the  use  of  the  Operational-J  parameter. 

3.4  THE  METHODOLOGY,  THE  APPROACHES,  AND  SOME  ISSUES 

In  the  previous  subsection  the  known  attempts  made  to  extend 

fracture-mechanics  representation  of  FCGR  in  the  elastic-plastic 

regime  were  presented.  Correlation  with  the  Operational  AJ 

2S  71—74  78  79  81 

parameter  seems  to  be  the  most  popular  among  authors  '  '  ' '  * 

elastic-plastic  FCGR  methodology.  As  their  results 
suggest,  AJ  has  correlated  FCGR  in  the  elastic-plastic  regime 
for  specimens  with  different  geometries.  The  correlation  is  very 
good  and  both  linear  elastic  and  elastic-plastic  data  lie  in  the 
same  scatter  band. 

One  can  argue  the  validity  of  using  AJ  as  the  correlating 
parameter  to  characterize  the  FCGR  in  the  elastic-plastic  regime, 
since  path  independence  of  the  J-tntegral  is  derived  for  a 
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(COD)-  (mm2) 


Figure  27.  High-Strain  Fatigue  Crack  propagation  *ates  as  a 
Function  of  Crack  Opening  Displacement,  Teats 
Carried  out  at  20°C  and  300°C  in  Repeated  Tension 
on  Mild  Stesel  Plates  1.83m  x  0.91m  x  12.7mm  Containing 
a  Central  Sharp  0*23m  slit  (Reference  82). 
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nonlinear  elastic  material.  When  the  material  undergoes 
unloading,  the  mathematical  validity  of  the  J-Integral  as  a 
parameter  which  characterizes  the  nonlinear  elastic  stresses  at 
the  crack  tip  is  violated. 

At  the  high  end  of  the  FCGR  curve,  the  crack  tip  moves 

forward  on  each  cycle  into  a  relatively  undamaged,  virgin 

material  in  terms  of  plastic  deformation  compared  to  the  intense 

deformation  zone  at  the  crack  tip.  Hence,  during  the  next  cycle 

for  which  AJ  is  calculated,  past  history  may  not  be  significant 

.  13 

compared  to  loading  .  It  must  be  noted  that  the  J-Integral  is 
derived  for  a  stationary  crack.  Figure  21  describes  two  different 
ways  of  defining  AJ  for  studying  FCGR  behavior.  To  define  AJ  for 
an  actual  component,  it  may  require  the  extensive  use  of  numerical 
analysis  and  possibly  the  contour  integral  definition  of  J. 

3.4.1  Life  Estimation  Scheme  Using  a  EPFM  Approach 

6  2 

Wilhem  et  al.  recently  presented  a  modified  stress 
intensity  factor  approach  based  on  elastic-plastic  J-Integral 
calculations  for  predicting  constant  amplitude  crack  growth  lives 
for  cracks  growing  through  nonlinear  stress  fields  established  by 
stress  concentrations.  Two  heat  treatment  levels,  -T651  and  -T7351, 
for  the  7075  aluminum  alloy  were  considered.  To  calculate  crack 
growth  lives  for  thru-thickness  cracks,  the  following  steps  are  taken 

1.  Perform  a  nonlinear  analysis  of  the  structure 
under  consideration  using  finite  element  analysis,  assuming 
Prandtl-Reuss  material  behavior. 

2.  Obtain  J  „  values  for  various  applied  stresses 

max 

and  crack  lengths. 

3.  Compute  modified  stress  intensity  factors  using 


KImax 
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plane  strain 
plane  stress 


(16) 
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Kl 


The  normalized  stress  intensity  factor 


max 


amax/rTa 

is  not  merely  a  function  of  crack  length  for  a  given  geometry 
as  in  LEFM  behavior  but  also  depends  on  the  remotely  applied 
stress  above  a  given  stress  level. 

4.  Using  KImax  computed  via  Equation  16,  a  typical  LEFM 
fatigue  crack  growth  rate  description  such  as  Fitzgerald  proposes 


da 

dN 


Km  (K  +  K  )  ( 1-R)  2 

max  '  max  env 


(17) 


can  be  used  to  calculate  FCGR  for  a  given  applied  stress  cycle. 


For  the  part-thru-crack  (PTC)  geometries,  the  authors  make 
the  following  assumption:  the  nonlinear  effects  (or  enhancement) 
for  both  thru-thickness  and  PTC  geometries  will  be  the  same  at 
similar  values  of  the  computed  elastic  stress  intensity  factor 
(K  ) .  For  this  reason,  it  is  suggested  that  the  (plasticity) 
modified  stress  intensity  factor  (K  )  be  written  as: 

r 


K 


P 


K 


4>  (a,  Ke)=  Ke 


(18) 


The  value  of  <f>  (a, Kfi)  is  a  correction  to  be  applied  to  the  elastic 
stress  intensity  factors  that  accounts  for  the  influence  of 
plasticity;  it  is  a  function  of  applied  stress  (a) ,  crack  length, 
and  the  plasticity  properties  of  material,  but  is  assumed 
independent  of  geometry.  It  is  not  certain  at  this  time  whether 
such  an  approach  will  provide  better  cross  correlations  than  those 
developed  using  the  LEFM  approach  modified  to  account  for 
plasticity  via  an  effective  crack  length  (real  length  plus  plastic 
zone  size) . 


3.4.2  Summary 


Most  of  the  available  FCG  data  generated  under  the 
elastic-plastic  conditions  are  presented  as: 

=  f  (Aj) 

A  limited  number  of  hysteresis  loop  diagrams  and  crack  length 
versus  number  of  cycle  curves  are  given  in  the  available  litera¬ 
ture.  Out  of  this  limited  information,  if  numerical  methodology 
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could  be  developed  to  calculate  other  possible  parameters ,  one 
could  test  the  available  data  with  these  parameters. 

5 

The  data  which  are  available  in  Dubensky's  paper  are  given 
as 


*  f  <AK)  (20) 

and 

a  =  g  (N)  (21) 

In  this  case,  as  there  are  no  hysteresis  loops  available,  It  may 
be  harder  to  develop  numerical  techniques  tc  evlauate  the  ability 
of  other  possible  parameters  which  could  correlate  PCGR  data.  It 
is,  therefore,  suggested  that  future  investigations  in  the  high 
FCGR  regime  report  displacement  parameter  information  as  a  function 
of  both  cycle  and  crack  length  histories,  as  well  as  the  raw 
cycle,  crack  length  data. 


SECTION  IV 


RECOMMENDATIONS 

Based  on  the  results  cf  the  Literature  review  presented  in 
Sections  1  through  3,  the  following  recommendations  for 
additional  work  are  made: 

1.  Conduct  additional  FCG  data  collection  studies  to 
determine  the  specific  set  of  conditions  -  material,  geometrical, 
and  mechanical  -  for  which  stress  and  structural  geometry  correla¬ 
tions  of  FCGR  by  the  LEFM  parameters  fail.  Two  separate  types  of 
studies  will  be  necessary: 

(a)  Study  FCGR  behavior  for  crocks  at  high  net 
section  stresses 

(b)  Study  FCGR  behavior  for  small  cracks  imbedded  in 
a  plastic  zone  at  the  root  of  a  notch. 

2.  Define  a  set  of  four  EPFM  parameters  for  further 
indepth  studies  to  evaluate  their  ability  for  correlating  FCGR 
data.  Delineate  the  advantages  and  disadvantages  of  each 
parameter.  The  set  of  EPFM  parameters  studied  should  include: 

(a)  A  LEFM  derivative  (e.g.,  K£  =  eB/ifa) 

(b)  The  Operational  J  (AJ  =  ~  ) 

B  3a  6 

(c)  The  Crack  Tip  Opening  Displacement,  <$t 

(d)  A  process  zone  parameter  (e.g.,  ) 
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